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^ Abstract 

Initially this paper reviews the mathematical issues surrounding the hydrostatic (HPE) and 
non-hydrostatic (NPE) primitive equations that have been used extensively in numerical weather 
prediction and climate modelling. Cao and Titi (2005, 2007) have provided a new impetus to this 
by proving existence and uniqueness of solutions of viscous HPE on a cylinder with Neumann-like 



this paper is whether extreme events are allowed to arise spontaneously in their solutions. Such 
events could include, for example, the sudden appearance and disappearance of locally intense 
fronts that do not involve deep convection. Analytical methods are used to show that for viscous 
HPE, the creation of small-scale structures is allowed locally in space and time at sizes that scale 
inversely with the Reynolds number. 

A review dedicated to the memory of Robin Bullough. 
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1 Review of the hydrostatic & non-hydrostatic primitive equations 

The hydrostatic primitive equations (HPE) were developed more than eighty years ago by Richardson 
as a model for large-scale oceanographic and atmospheric dynamics [T]. In various forms they have 
been the foundation of most numerical weather, climate and global ocean circulation predictions for 
many decades. The HPE govern incompressible, rotating, stratified fluid flows that are in hydrostatic 
balance. This balance is broken, however, in deep convection which occurs in cloud formation, flows 



over mountains and vertical fluid entrainment by strong gravity currentsM In this context, the 



hydrostatic approximation is the most reliable and accurate of all the assumptions made in applying 
them to operational numerical weather prediction (NWP), climate modelling and global numerical 
simulations of ocean circulation. The approximation arises from a scale analysis of synoptic systems; 
see I21[31|4l[5l[a[71[ai3[Ta[TTl[Tl[T31[I3[I3[ia[T3 Basically, it neglects the vertical 
acceleration term in the vertical momentum equation. Physically, this enforces a perfect vertical 
force balance between gravity and the vertical pressure gradient force, which means that the pressure 
at a given location is given by the weight of the fluid above. 

Until recently, the resolution of operational NWP models has been limited by computing power 
and operational time constraints on resolutions in which the hydrostatic approximation is almost 
perfectly valid. Operational forecasts therefore have relied mainly on hydrostatic models, which 
applied extremely well in numerical simulations of the global circulation of both the atmosphere 
and the ocean at the spatial and temporal resolutions that were available at the time. In parallel, 
the development of non-hydrostatic atmospheric models, which retain the vertical acceleration term 
and thus capture strong vertical convection, has also been pursued over a period of more than 
four decades, especially for mesoscale investigations of sudden storms. As computer systems have 
become faster and memory has become more affordable over the past decade, there has been a 
corresponding increase in the spatial resolution of both NWP and climate simulation models. This 
improvement has facilitated a transition from highly developed hydrostatic models towards non- 
hydrostatic models. Over the past decade, atmospheric research institutions have begun replacing 
their operational hydrostatic models with non-hydrostatic versions [I4l [T5l [17] . 

An impending 'model upgrade' may be tested by comparing the solution properties of the existing 
hydrostatic model at a new higher resolution with those of its non-hydrostatic alternative (NPE). Us- 
ing the fastest supercomputers available, global simulations of atmospheric circulation at resolutions 
of about 10 km in the horizontal (in the region of the hydrostatic limit) have recently been performed 
[19l[2Q]. In the past these simulations have been carried out using only hydrostatic models, as the 
non-hydrostatic versions were still under development. Perhaps not unexpectedly, the simulations at 
finer resolution found much more fine-scale structure than had been seen previously. More impor- 
tantly, the smallest coherent features found at the previous coarser limits of resolution were no longer 
present at new finer resolutions, because the balance between nonlinearity and dissipation that had 
created them previously was no longer being enforced. Instead, it was being enforced at the new 
limits of resolution. 

^The subtle difference is explained in fj2]between what is known as modellers' HPE (MHPE) pioneered by Richardson 
and incorporated in many numerical weather prediction and climate simulation models, and HPE in cartesian geometry 
discussed here. 
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Consequently, an important conclusion from these PE simulations has been that the transition 
from the existing hydrostatic model to its non-hydrostatic alternative, and the distinction between 
their solutions, will require in both cases a much better understanding of the formation of fine-scale 
frontal structures [221121]. For discussions of the history of ideas in fronto-genesis see [221 123 1251 1211 ■ 
An example of fine-scale wind measurements is a 'gust front' are found in Figure 1 [TBI [2S1 [2Z]- Gust 
fronts are indicated by 'bow echoes' on Doppler radar. These have been extensively observed and 
studied as the precursors of severe local weather. 

In particular, the physical parameterizations in hydrostatic models are likely to require dramatic 
changes when simulated at new finer resolutions. Therefore, in order to obtain the full benefit of the 
potential increase in accuracy provided by such high spatial resolution when using the non-hydrostatic 
model, one must first determine the range of solution behavior arising at the finer resolution in the 
computations of the previous hydrostatic model [20|. Over relatively small domains, operational 
model resolutions in the atmosphere are in fact already beyond the hydrostatic limit. These mod- 
els run at resolutions finer than 10 km, where convection is partially resolved (thus, inaccurately 
simulated) by non-hydrostatic models. However, convection can only be fully resolved at model res- 
olutions of about two orders of magnitude finer, that is, about 100 m or less in the horizontal. It is 
therefore likely that, for many years to come, operational hydrostatic and non-hydrostatic models will 
function at resolutions where convection can only be partially resolved. Consequently, the distinction 
between hydrostatic and non-hydrostatic solution behavior at these intermediate resolutions becomes 
paramount for the accuracy, predictability, reliability and physical interpretation of results in NWP. 



Figure 1: An example of a gust front manifesting itself in a sudden rise in wind speed [161 126| . 

Against this background of renewed interest in the improvement of hydrostatic models and their 
newly acquired computability at higher resolution, this paper also reviews the development of new 
mathematical ideas introduced by Cao and Titi |L28ji29j, who have proved the existence and unique- 
ness strong solutions of the viscous hydrostatic primitive equation^ The corresponding regularity 
problem for viscous NPE remains operrj Ju [33J then built on the result of Cao and Titi to show 




^See also Kobelkov j30l 1311 [32] for a subsequent independent method. 

^This is closely related to the open regularity problem for the three-dimensional Navier- Stokes equations which has 
been solved only for very thin domains (0 < Qa 1) : see [ 341 1351 [36] . 
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that the viscous HPE actually possess a global attractor, which means that the bounds are constants 
although no explicit estimates were given. The status of results prior to the Cao and Titi proof 
[221 [22] (e.g. the existence of weak solutions) can be found in the papers by Lions, Temam and 
Wang [371 ES ESI and Lewandowski [4D]. Considering the technical difficulties of this new work, 
which has been solved on a cylindrical domain with Neumann-type boundary conditions on the top 
and bottom, these results must be classed as a major advanc^ The significance of this result in 
establishing the existence and uniqueness of strong solutions of viscous HPE is that no singularities 
of any type can form in the solutions. 

The title and abstract of this review refer to the idea of 'extreme events' : this term should 
not be interpreted in the present context in a statistical sense, but as the spontaneous appearance 
and disappearance of local, sudden, intermittent events with large gradients in the atmosphere or 
the oceans. These events would register as spatially localized, intense concentrations of vorticity 
and strain. Until the the Cao-Titi result, a major open question in the subject had been whether 
these strong variations ultimately remained smooth? Now that we know that solutions of HPE 
must remain smooth at every scale, the question still remains whether this smoothness precludes the 
existence of extreme events representing front-like structures: this would be manifest in jumps or 
steep gradients in double mixed derivatives such as \d'^ui/ dx2dx^\ or \d'^U2/dx\dxz\ within finite 
regions of space-time. The ultimate aim is to explain fine-structure processes such as those that 
occur in fronto-genesis : see the work of Hoskins [23 121] and Hoskins and Bretherton [25] for 
general background on this phenomenon. 

f|3] of this paper seeks to address this problem from a mathematical angle. The regularity of 
solutions opens a way of showing that front-like events are possible at very fine scales. The method 
used is to search for fine structure by examining the local behaviour of solutions in space-time. The 
conventional approach to PDEs has been to prove bounded-ness of norms which, as volume integrals, 
tend to obscure the fine structure of a solution. The idea was first used (under strong regularity 
assumptions) to estimate intermittency in Navier-Stokes turbulence [41]. The mathematical approach 
is to show that solutions of HPE may be divided into two space-time regions S+ and S^. If is 
non-empty then very large lower bounds on double mixed derivatives of components of the velocity 
field (ui, U2, u^) such as \d'^ui/dx2dx^\ or \d^U2/dxidx3\ may occur within S~. These large point- 
wise lower bounds on second-derivative quantities represent intense accumulation in the {xi,xs)- and 
(x2, a;3)-planes, respectively. Clearly, if solutions of great intensity were to accumulate in regions of 
the space-volume, this could only be allowed for finite times, in local spatial regions. Thus one would 
see the spontaneous formation of a front-like object localised in space that would only exist for a finite 
time. More specifically, the lower bounds within referred to above are a linear sum of 7^^^^ , TZ^^ 
and TZ^ j< where Tiuh^r ^'^^ IZu^ are local Reynolds numbers. That is, they are defined using the local 
space-time values of Uhor{xi,X2, x^, r) = \/u\ + and the vertical velocity = u^{xi,X2,x^^ r). 
Likewise, TZa,T is a local Rayleigh number, defined using the local temperature T{xi, X2, X3, r). 
In contrast. Re is the global Reynolds number. These lower bounds within the S^-regions can be 
converted into lower bounds on the two point-wise inverse length scales \h 

LX^ > Cu^^^TZuhor + CuaTZua + Ci,T'^a,T + C2,T-Re^/^7^^^^ + forcing . (1.1) 

* Unfortunately the methods used for HPE do not extend to the three-dimensional Navier-Stokes equations; although 
HPE had been thought to be the more difficult of the two. See Lions, Temam and Wang [37] l38l |39] . 
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The term TZu^ ~ 0{e) is negligible compared to TZuhor because regions of strong vertical convection 
are absent in HPE. If such a result were valid for NPE, it would be this vertical term that would be 



restored. The estimate for Xh in (1.1) is of the order of a metre or less, at high Reynolds numbers. 
Our interpretation of this very small length is not that it is necessarily the thickness of a front, but 
that it may refer to the smallest scale of features within a front. If an equivalent result for NPE 
could be found then the estimate of the length-scale Ajv (say) would likely be considerably smaller 
than that for A^- 



2 Formulating the Primitive Equations 



Quantity 


Symbol 


Definition 


Typical horizontal length 


L 




Typical vertical length 


H 




Typical temperature 


To 




Aspect ratio 


OLa 


aa = H/L 


Typical velocity 






Twice vertical rotation rate 


f 




Global Reynolds number 


Re 


Re = UoLv-^ 


Rossby number 


e 


e = UoiLf)-^ 


Rayleigh number 


Ra 


Ra = gaToH^uK)-^ 


Hydrostatic const 


ao 


ao = eaa^'^RaRe^^ 


Frequency 


Uo 


u)o = UoL-^ 



Table 1: Definition of symbols for the primitive equations. 



In what follows, dimensionless co-ordinates are denoted by {x, y, z, t). These are related to di- 
mensional variables {xi, X2, x^,t) through the horizontal length scale L, the vertical length scale H : 
Table 1 gives the various standard definitions. Table 2 gives the notation concerning dimensionless 
and dimensional variables and some of the relations between them. 

Dimensionless versions of both the HPE and NPE are expressed in terms of a set of velocity 
vectors based on the two horizontal velocities (u, v) and the vertical velocity w. These form the 
basis for the three-dimensional vectoi0 

V{x,y,z,t) = {u, V, ew) , (2.1) 

satisfying div V = or Ux + Vy = —ewz- The hydrostatic velocity vector is defined by 

v = {u,v,Q), (2.2) 

and the non-hydrostatic velocity vector by 

= (u, V, a^a^w) , (2.3) 
^Comparing with the notation in |9j, {u, v) = 1x2 while V = U3 and D = V3. 
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Dimensionless 


Dimensional 


Relation 


Horiz co-ords 


x,y 


Xi,X2 


(x, y) = {xi,X2)L-^ 


Vertical co-ord 


z 


X3 


z = X3H-^ 


Time 


t 


T 


t = tUqL-^ 


Horiz velocities 


{u, v) 


u = (ni, U2) 


(n, v) = {ui,U2)Uq^ 


Vertical velocity 


w 


U3 


eUoaaW = U3 


Velocity vector 


V = {u, V, ew) 






Hydro- Velocity 


V = {u, V, 0) 






Non-Hydro-Vel 


= {u, V, a'^ew) 






Temperature 


e 


T 


e = tTq^ 


2Z)-Gradient 


V2 = id^, + jdy 






SD-Gradient 


V3 = id:^ +jdy + kdz 


V = id^^ + jd^2 + kaadxs 


V3 = LV 


3D-Laplacian 


A3 = a2 + a2 + 52 






Vorticity 


(jL> = curl V 


ft 


ila = UJ3UJ0 


Hydro-vorticity 


C = curl V 






Non-hydro-vorticity 


to = curl 






Heat transport 


q 


Q 





Table 2: Connection between dimensionless &i dimensional variables. 



where aa is the aspect ratio defined in Table [3] Neither v nor d are divergence-free. 

The HPE that are studied in this paper should be compared with the modellers' HPE (MHPE) 
used by Richardson and incorporated in many numerical weather prediction and climate simulation 
models to this day [HI [151 IISI [U]- While the MHPE are written for a compressible fluid in the 
non-Euclidean geometry of a shallow atmosphere (i.e. the domain of flow is a spherical shell of 



constant radius, although vertical extent and motion are allowed) the HPE as laid out in Section 2.1 
are written for an incompressible fluid in a Cartesian geometry. The 'incompressibility' of hydrostatic 
models appears when pressure is used as vertical coordinate. Hoskins and Bretherton [25] used a 
function of pressure to define a vertical coordinate that delivers the hydrostatic primitive equations 
for a compressible fluid that are nearly isomorphic to those for an incompressible fluid when ordinary 
height is used as vertical coordinate. They then introduced an approximation in the continuity 
equation that made the near isomorphism exact. The resulting HPE have frequently been used in 
theoretical studies and have led to many important results, but they are not identical to the MHPE 
that numerical modellers since Richardson have typically used. 

2.1 Hydrostatic primitive equations (HPE) 

The HPE in dimensionless form for the two horizontal velocities u and v are 

e(^^^ + V ■V,}ju-v = eRe-^A3U-p^, (2.4) 



30th/04/10 J D Gibbon and D D Holm 



Extreme events in climate models 



7 



and 

^^ + V ■V3^v + u = eRe-'Asv-py. (2.5) 

There is no evolution equation for the vertical velocity w which lies in both V • V3 and the in- 
compressibility condition d'wV = 0. The z-derivative of the pressure field p and the dimensionless 
temperature B enter the problem through the hydrostatic equation 

ao@ + Pz = 0, (2.6) 



where ao is defined in Table 1. The hydrostatic velocity field v = (n, v, 0) appears when (2.4), (2.5) 



and (2.6) are combined 

'd 



Q^+V ■ V3 j v + k XV + aoke = ei^e"^ A3 v - Vsp , (2.7) 

and taken in tandem with the incompressibility condition divt; = —ewz- The dimensionless temper- 
atur^©, with a specified heat transport term q{x,y, z,t), satisfies 



+ F • V3 j e = {aRey^A^e + q . 



(2. 



The chosen domain is a cylinder of height H {0 < z < H) and radius L, designated as C{L,H). In 
terms of dimensionless variables, boundary conditions are taken to be: 

1. On the top and bottom of the cylinder z = and z = H the normal derivatives satisfy 
Uz = Vz = Qz = together with w = 0. This is appropriate for the atmosphere where the 



heat transport term q{x, y, t) in (2.8) is dominant. Ideally, in the case of the ocean, as in 
[211 ES]. the zero conditions for and Vz on the cylinder top would need to be replaced by a 
specified wind stress, but this creates surface integrals which can only be estimated in terms 
of higher derivatives within the cylinder. 

2. Periodic boundary conditions on all variables are taken on the side of the cylinder Sc, which 
is a minor change from those in [2S1I2S1E3]. 

The key feature for HPE is the use of enstrophy (i.e. vorticity squared) rather than kinetic energy as 
a quadratic measure of motion. This means that the pressure gradient is excluded from the dynamics 
at an early stage and, more importantly, that the analysis is immediately in terms of spatial derivatives 
of velocity components rather than the components themselves. In contrast, earlier approaches based 
available potential energy include the pressure field until a global spatial integral removes it and this 
has made local results more difficult to reach. In addition, energy arguments deliver results about 
attainable velocities not about velocity gradients. 
The vorticity C = cur\v is specifically given by 

^ = curl V = -ivz + juz + k{vx - Uy) . (2.9) 



®Note that we use the upper-case O for the temperature to avoid confusion with lower-case 9 which is conventionally 
used for the potential temperature. 
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This contains the same fc-component as the full 3D-vorticity cj = V3 x V but with the i«-terms 
missing. The relation between ^ and u gives rise to a technically important relation 



Taking the curl of (2.7) gives 



dC 

e-^ = eRe'^AsC + ecur\{V X C) - cur\{k X V + aoke) , 
This relation forms the basis of the proof of Theorem [1] in ^JS) 

2.2 Non-hydrostatic primitive equations (NPE) 



(2.10) 



(2.11) 



The NPE restore vertical acceleration so that, in contrast to (2.6), the equation for w reads 



d 



ale^ [j^^+V ■V:i]w + aoQ + p^ = eRe 'A3W . 



Using the definition of t) in (2.3) and putting (2.12) together with (2.4) and (2.5) we find 

e + V ■ V3 j + fc X D + aoA;9 = ei?e"^ A3 1) - Vsp . 



(2.12) 



(2.13) 



The equation for the temperature (2.8) remains the same. In contrast to the hydrostatic case where 



regularity has been established with realistic Neumann-type boundary conditions on the top of a 
cylinder, several major open questions remain in the non-hydrostatic case : 



1. No proof yet exists for the regularity of the NPE (2.13). Their prognostic equation for w brings 



NPE much closer to the three-dimensional Navier-Stokes equations for a rotating stratified 
fluid than HPE, in which w is diagnosed from the incompressibility condition. While regularity 
properties for Navier-Stokes fluids on very thin domains (0 < aa « 1) has been established 
[Ml [35], the general regularity problem remains open. 



2. The boundary integrals for HPE estimated on the top and bottom of the cylinder C{H,L) are 
also problematic for NPE. Whereas the vorticity in the hydrostatic case has only terms Uz 
and Vz in the horizontal components, which are specified on the top and bottom of C{H,L), 
the non-hydrostatic vorticity 

to = curl t) 



(2.14) 



also includes Wx and Wy terms which are not specified. 



Similar to (2.10), there exists a technical relation between V, and tv 

- F . V30 = V X tt) - iV3 [n2 - eW{al - if] , 



so the curl-operation on (2.13) gives 
dtv 



dt 



eRe ^A^to + ecurl [V x to) - curl {k x ti + aokQ) 



(2.15) 



(2.16) 



30th/04/10 J D Gibbon and D D Holm 



Extreme events in climate models 



9 



without involving pressure terms. 



To address the influence of the material derivative in (2.12), typical values of a^e can be 



determined. Typical observed values for mid-latitude synoptic weather and climate systems are: 



Ota 

w/u 

£ 



H/L ^ lO^m/lO^m 



10 



-2 



10 



-3 



(2.17) 



U/{foL) ^ Wms^^/iW^s-'^W^m) « 10"^ 
For mid-latitude large-scale ocean circulation, the corresponding numbers are: 



"a 

W/U 

£ 



H/L lO^m/loV ^ 10 



-2 



,-1 



10 



-2 



lO^'^ms^VlO^^ms 
U/{foL) ^ IQ-^ms-'^ /{IQ-'^s-^lO^m) ^ 10" 



(2.18) 



Thus, for these typical conditions, a^e^ ~ 10"*^ — 10~^ <C 1 so the hydrostatic approximation can be 
expected to be extremely accurate in calculations of either synopic weather and climate, or large-scale 
ocean circulation at mid-latitude: see [4511431144]. 

3 Extreme events 8^ their interpretation 

The proof of global existence and uniqueness of solutions of HPE by Cao and Titi guarantees that at 
each point in space time a unique solution exists (2S1I2S1E3]- '^^^^ result can be exploited to prove 
Theorem 1 and from this to consider the idea of extreme events. 



HPE quantity 


Definition 


Local (horizontal) Reynolds no in terms of \u\ 


^hor 1 1 


Local (verterical) Reynolds no in terms of \us\ 


IZu.^ = L Ms z^^i 


Local Rayleigh no in terms of T 




/3i 




/?2 






o 9 I 243 


Pus 


^"3 = ( + ll^l) 


Pl,T 


/3l,T = f (1 + ^) ^a'^f 


1^2, T 




V, 





Table 3: Definitions in the hydrostatic case. 



Consider the arbitrary positive constants < ^2 < 1 chosen such that 

253 = l-|<5i-52 >0. 



(3.1) 
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1 

252 



"'ho 



+ 



a 



Re\Q(' 



(3.2) 



These are used to define a hydrostatic forcing function Fh[Q) 

Fh{Q) = e-^ {52 

The spatially global quantity 'H(i) 

^(t)= / {Re^\C\^ + Re\e^\^)dV . (3.3) 
Jv 

is now used in the proof of the following theorem : 

Theorem 1 In C{H, L) over some chosen time interval [0, r*], the space-time 4-integral satisfies 



-Pi 



+ 



|V2J^3p + a« 





2 

+ «a 


8'^U2 


2' 


8x28x3 


8x18x3 





l32al\V3T,f + (3^,K^^ 



(3.4) 



where the coefficients are given in Table\^ 



Remark 1 : The proof in <:3.1 shows that the right hand side of (3.4) is, in fact, \%{t^,). Because 
solutions of HPE are regular this is bounded above for all values of t^. However, it has a lower bound 
of zero which, although not necessarily a good lower bound, is uniform in r*. Regularity of solutions 



is also a necessity in order to extract point-wise functions from within the 4-integral in (3.4). 

Remark 2 : Of course, any integral of the form ff{A — B) dVdt > trivially indicates that there 
must be regions of space-time where A > B but, potentially, there could also be regions where 
A< B. Such integrals are common, particularly energy integrals, and in most cases the information 



gained is of little interest. In the case of (3.4), however, the proof in 3.1 will demonstrate that much 



effort has gone into rigorously manipulating it into a form that produces sensible and recognizable 



physics. From (3.4) we conclude that: 



1. There are regions of space-time S+ C on which 



/3i 



|V2^^3|^ + a 



8'^ui 


2 

+ al 


d^U2 


T 


8x28x3 


8x18x3 





p2al\VT^., 



< Pu,^nt,^^ + Pu,K^ + Pi,t111t + P2,TRe^nlT + Fh{Q) + O(r-i) (3.5) 
2. Potentially there are also regions of space-time C on which 



/3i 



8'^ui 



8x28x3 



d'^U2 



8x18x3 



+ (32al\VT,,\ 



> (^u,^rK.or + f^n.K, + /3l,T7^^,T + P2,TRe'nlT + Fh{Q) + ©(t'^) (3.6) 
It must be stressed that there is no information here about the nature of the two sets : the 



4-integral in (3.4) above yields no more information other than the possibility of a non-empty set 
existing. It says nothing about the spatial or temporal statistics of the subsets of (which 
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may have a very sensitive r*-dependence) nor does it give any indication of their topology. These 
results, however, are consistent with the observations of fronts in the atmosphere where 
large second gradients appear spontaneously in confined spatial regions often disappearing 
again in an equally spontaneous manner. This behaviour would also be consistent with S~ being 



comprised of a disjoint union of subsets although no details can be deduced from (3.4) 



To illustrate the nature of a front, very large values of double mixed-derivatives are required in 
local parts of the flow as envisaged in the early and pioneering work of Hoskins |231 For instance, 
very large values of \d'^ui/ dx^dx2\'^ or \d'^U2/ dx'idxi\'^ would represent intense accumulation in the 
(xijXs)- and (x2, X3)-planes respectively. 

The large lower bounds in Theorems [T] can be interpreted in terms of a length scale. To achieve 
this, define the point-wise inverse length scale A^^ such thal[^ 



-n6 



r2, ,-2 



dx2dx3 



+ a„ 



d^U2 



dxidxs 



+all32(3^'L^\V3T,,\-'T^\ (3.7) 



The sixth-powers on the right hand sides of Theorem [T] gives results within S 

LXJj^ > Cu^^^TZuhor + CusTlus + ci^xT^a^T + C2,T^e^/^7?.y^ + forcing 



(3.8) 



where the coefficients c^j^^,,, c^g and q^t can be calculated from Theorems[T]: indeed the term TZ. 



can be ignored as the vertical velocity U3 ~ 0{e). (3.8) can be interpreted as a lower bound on the 



inverse length scale of the smallest feature in a front. As already noted, no information is available 
on any statistics nor on the shape and size of subsets of S~. 

3.1 Proof of Theorem [T] 

3.1.1 The evolution of the enstrophy 



The equation for the HPE enstrophy C in (2.11) is now considered on the cylinder C{H,L) with the 
boundary conditions Uz = Vz = w = on the top and bottom but with periodic side-wall conditions 



/ ICl'^dV = eRe-^ [ C'^sCdV + e [ C ■ cur\ {V x C) dV 
dt Jv Jv Jv 

- / C-curl(fc X i7-Faofce)dy. 
Jv 

The notation in the rest of the paper is 



\f\PdV 



i/p 



(3.9) 



(3.10) 



Note that in estimating the integrals in (3.9), a standard vector identity and the Divergence Theorem 



are used in which a surface integral naturally appears each time. However, on the cylinder top and 
bottom 

fcxC = on z = 0, H . (3.11) 



^The point-wise local length scale Ajj^ is formed in the same dimensional manner as the Kraichnan length £k, 



namely, from a combination of the palenstrophy jVajj and the viscosity u given by ^j. 
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This and the side-wall periodic boundary conditions make zero all the surface integrals that appear 
from the Divergence Theorem. 

Considering (3.9 term by term, the first is (div(^ = 0) 

/ C-A3CdV = - [ C -curl curl C dT/ = - / |curl f^l/" • (3.12) 
Jv Jv Jv 



The second integral on the right hand side of (3.9) 



/ C ■ cur\ {V X C)dV = [ curl C ■ {V X C)dV , (3.13) 
Jv Jv 

and so 

[ C-cur\{V xC)dV < llcurlClbllVllellClls. (3.14) 
Jv 

Thirdly, 

C ■ cur\{k X V + aoke) dV = cur\ C ■ {k x v + aoke) dV . (3.15) 
Jv Jv 

Thus (3.9) can be re-written as 

le^ [ ICl'^dV < -eRe-^ [ IcurlCp + ellcurl ClbllV^I 
dt Jv Jv 



I6||t,| 3 



- [ {curl C) ■ {k X V + aokQ) dV 
Jv 

= -eRe-^ [ IcurlCp + Tl + T2 + T3 , 
Jv 



(3.16) 



where 

Tl = e||curlC||2||V||6||C||3 (3.17) 

J2 = - [ (curIC) ■{kxv)dV , T3 = -qq [ (curIC) • (kQ) dV . 

Jv Jv 

Before estimating Tl, T2, and T3, we prove 



(3.18) 



Lemma 1 In the cylinder C{H, L) 

||C||3<3V2||cor/C||^/'||^||^/'. (3.19) 

Proof: Consider 

/ Kl^dV = [ C-(CC)d^ 
Jv Jv 

= [ v-cur\{CC)dV- [ d\y{CivxC))dS 
Jv Js 

= [ V (Ccurl C + (VaC) x C)dV - [ div {({v x C))dS . (3.20) 
Jv Js 

The surface integral is zero because fc x = on the cylinder top and bottom. A vector identity 

VsC = |(C)~'V3(C • C) = C • VaC + C X curl C (3.21) 
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allows us to write 



curl (CO = Ccurl C + (C x curl C) x C - C x (C • V3C) , 



so we conclude that 



\C\^dV= / V- CcurlC + (C X curIC) x C-C x (C • V3C) dV . 



Using a Holder inequality, it is then found that 

/ |C|=^tZT^<3||curlC||2||CI|3||^^| 

□ 



giving the result (3.19) 



(3.22) 



(3.23) 



(3.24) 



Lemma 2 Within the cylinder C{L, H), Tl, T2 and T3 are estimated as 

\T1\ < ss^eRe-'\\curlCf2 + ^^Re^2\\V\\l + \\v\\l) 

for any di > 0. Moreover, for any ^2 > and 82 > Q 

\T2\ < i(52ei?e~icur/C||i + i<5^^e"^i?e||i;||i, 



I 73| < i62eRe-'\\curl CWl + ^^M&Wl ■ 



(3.25) 

(3.26) 
(3.27) 



Proof : In the following the 6i > are constants introduced by a series of Young's inequalities; 
1) Using Lemma [1]T1 can be written as 



Tl = ellcurlClbll^l 



1/2,, ,,1/2 



< e||curlC||2||V||6 X 3i/2||curl Clla^^H-!; 

< (<5iei?e-icurlC||i)'/'(9e5r'^e3||F||6)i/^(9e5r=^^e3||^||6)^/^2 

< f <5ieiie-i||curl C||i + ^^Re' {2\\V\\I + \\v\\l) . 



2) From (3.18), T2 can be estimated as 

|T2| < llcurlClbll-J^lb = {52eRe~^\\cur\C\\2)^^^ {62^e-^Re\\v\\l)^^^ 
< i(52ei?e~^||curlC||2 + ^S^^e-'^ReWvg. 



(3.28) 



(3.29) 



(3.30) 



3) From (3.18), and using the same constant 82, T3 is estimated as 



|T3| < ao||curlC||2||e||2 = (52ei?e-i||curlC||i)'/'(^^i?e||G||i 
< |<52ei?e-i||curlC||^ + ^i?e||G||l. 



1/2 



(3.31) 



as advertized. 



□ 
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Returning to (3.9), a division by e and a gathering terms gives 



\jJ^\C?dV < -{l-l6i-62)Re-' j^\cuv\C\^dV (3.32) 
+ Ai^e^ (2e6||^||6 + 3||^||6) + ^(^62 + Re\\v\\l + 2^^e||e|| 



2 

2 • 



The following Lemma relates \cur\<^\^dV to sums of squares. 

Lemma 3 Let uj-^ be the third component of the full vorticity cj = V3 x V = (uji, uj2, wa). Then 
for any < 60 < 1 



[ \curlC\^dV > min{l, 2(l-(^o)} / \V 2i03\^ dV + 26o [ {\uy,\^ + \v,,\^}dV 
Jv Jv Jv 

+ / {\uzz\'^ + \v,^\'^ + e^{l-6o)\w,,\'^]dV . 
Jv 



(3.33) 



Proof : Using C, in (2.9) in the form C, = —Vzi + Uzj + ujsk with = Vx — Uy and recalling that 

ewz = -{ux + Vy) 



curl C, 



i j k 

dx dy dz 

-Vz Uz UJ3 



i{uJ3,y - Uzz) - j{i^3,x + Vzz) " ekWz 



Thus 



Icurl Cr = |V2t^3| + {Uzz + Vzz + e W^L) - 2(^3,5,^^^^ - W3,xf 



(3.34) 



(3.35) 



Now, invoking the boundary conditions and the fact that L03 = Vx — Uy, integration by parts gives 



- 2 / {ujs^yUzz - ^3,xVzz)dV = 2 \uJs,z\'^dV 
Jv Jv 



and so 



/ Icurl C|'<iy= / {\V2^3f + 2\uj3,z\^dV}+ [ {ul + vl + eW,,}dV. 
Jv Jv Jv 

Moreover, it is easily shown that 

+ 2\u3,zf)dV = [ {{uxz-Vyzf + 2{ul, + vl,)]dV 

Jv 



(3.36) 



(3.37) 



I 2| |2 



> 2 / {uL + <z)dV 



(3.38) 



where a pair of horizontal integrations by parts in the first line of (3.38) have been performed. A 
linear combination of (3.37) and ( |3.38[ ) gives 

/ Icurl Ci'dy = / {\V2^3? + 2\u3,z?dV]+ [ {ul, + + sW,,}dV 
Jv Jv Jv 

> I {|V2C^3|' + 2(l-(5o)|a;3,.|2dF} 
Jv 

+ ^ {ul + vl + (1 - 6o)eW,,}dV + 26o j^iu^ + vl) dV , (3.39) 



which gives (3.33). This completes the proof. 



□ 



30th/04/10 J D Gibbon and D D Holm Extreme events in climate models 15 

3.1.2 The evolution of L \Qz\^dV 



The partial differential equation for given in (2.8) with BCs applied on C{L,H) 

+ V-Vse = {aRe)-^A3e + q (3.40) 

ot 

is now differentiated with respect to z to give 

= -{aRe)-^ I iVse.l'dF- / V-Vsi^Q^dV 
Jv Jv 

- [ e,iu,Q^ + v,ey + ew,Q,)dV + [ e,q,dV. (3.41) 
Jv Jv 

However, given that divV^ = and w = on 

fv-V3{\Ql)dV = [ {d\y{ielV)-IQld\yV}dV 
Jv Jv 

= 1 / {n.v)elds 

Js 

= ±|e / wQldxdy = 0. (3.42) 

Js± 



Integrating by parts the 3rd and 4th terms in (3.41) gives 



-{al 



+ / @{u^@xz + VzQyz + eWzQzz)dV - / QzzQdV 

Jv Jv 

TRe)-^ [ \V3ez\^dV+ [ e{uzQ,z + Vzeyz + ewzezz)dV 
Jv Jv 

QzzqdV (3.43) 
where div V = has been used. Using a Holder inequality it is found that 

\jj \Qz\''dV < -{aRe)-' j iVsQzl^dV+WQzzhhh 

+ l|0||6{||M.||3||e,,||2 + ||w.||3||ej,,||2 + £|h,||3||e..||2} . (3.44) 

The next task, addressed in the following Lemma, is to estimate ||uz||3, \\vz\\3 and \\wz\\3 in terms 
of their second derivatives. 

Lemma 4 With Neumann boundary conditions on C{H, L), the vector Vz and the scalar Wz satisfy 

\\vzh<Q^'^\\vzz\\T\\v\\]J\ (3.45) 
\\wzh<2"^\\wzz\\TM\T ■ (3-46) 
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Proof: 



[ \v,\^dV = [ {\u,\^ + \v,\^f/''dV 
Jv Jv 

\u,f + \v,f)dV 



(3.47) 



and, given the boundary conditions on u and v, 



\u,rdv 



U::Uz\Uz\dV 



UUzz\Uz \ + UU' 



d\uz 
dz 



dV 



uuz I Uz I dxdy 



(3.48) 



which holds because d\f\/dz < \fz\ for any appropriately function differentiable /. Thus (3.47) 
becomes 

/ \vz\^dV < 3 / {Inllii^^llii^l + 
Jv Jv 

< 6 / |t'||l»22||^'2|d^ 

Jv 



< G\\v\\e\\Vzz\\2\\Vz\\3 , 

which gives the advertised result. The result for w follows in a similar manner. 



(3.49) 

□ 



Continuing with (3.44), multiplying by Re'^ , where 7 is to be determined, (3.44) becomes 



I^Ee^ / IG.I^dF < -{aRe)-'Re^ [ \V3ez\^dV + 2^/hRe^\\eMw\\l/''\ 
dt Jv Jv 

+ 6i/2^e^||^,,||^/'||^||^/'||e||6{||e,,||2 + ||ej,,||2} 

+ i?e^||e,,||2|k||2. 
In turn, this re-arranges to 

ii?e^^ / \ez\'^dV < -{aRe)-^Re< I {VsQzl'^dV 

dt Jyi Jv 



W 



1 1/2 



Z2 112 



lie 



zz\\2 



(3.50) 



53Vi?e^M|0|l[ 



+ [463Re-'e^Wzz\\lY^^ [63{aRer'Re''\\ezz\\l 
1/6 r._, « .... „«li/i2 



21 1/2 



6^^e^a^Re^'\\w\ 



+ { [46,Re-'\\vzzf2Y^' [3%'Re^^v\\lY^''} 



161I/12 



1/6 



{26i{aRey^Re"' [\\e^z\\l + W&yzWl] } 



21 -11/2 



+ {S^iaRey'Re^WezzWlY^' {6^\aRe)Re^\\qg} 



12-1 1/2 



(3.51) 
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where 2bi + 62 = 9 + 67 and 2ci + C2 = 9 + 67. Using Young's inequality it is found that 

+ 5^Re-h^\\w,,\\l + \5^{aRe)'^Re^\\Q,,\\l 

+ \6fa^Re'-m\l + ^^67^h' a' Re'^w\\l 

+ 6iRe-^\\v,,\\l + 5A{aRe)-^Re^ [W^xzWl + \\Qyz\\^\ 



+ \8^{aRer^Re'\\Q,,\\l + \5l\aRe)Re^\\q\\l (3.52) 



Gathering terms we find 

\Re^j^ le.l^dF < -{aRe)-^Re< 5i){\Q^,\'' + \Qy,\'') + (1 - iJg - W)\Qzz?] dV 



12SI 



|6 



+ ii5-'(<Tfle)ife'|k||L (3,53) 

3.1.3 A combination of the fluid and temperature inequalities 
(3.53) is now combined with (3.32) 

< -{1- -^61-62) Re-^ [ \cur\C\'^dV 

Jv 

- {aRe)-'Re^ [ {{1 - 64){\e^z\^ + \Qy,\^) + (1 - iJg - |<55)|e..|2} dl^ 
Jv 

+ Re-^ {SaWvzzWI + he^w,,\\l) + ( Jg^e^ + ^«e^') ll^'lle 

+ H'^' ^ i) ^ ^ |55"Hai?e)iie^|k||i (3.54) 



and 



2bi + 62 

2ci + C2 



9 + 67 (3.55) 



Our choices are 

61 = ci = -3, 62 = C2 = 3, 7 = -2 . (3.56) 
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^3 and 64, also need to be chosen such that the ||i>z2|p- and e^Hty^^lP-terms cancel from ||curl(^||2. 
To this end we choose (5o = = and make 



{1 - I6i - 62) {1 - 60) = 63. 



(3.57) 



Hence 64, = 263. With this we choose S5 such that the coefficients 1 — 64 and 1 — |(53 — ^65 within 
the double derivatives of the temperature are equal. Thus ^5 = 3^3. Together we have 



'1 



263, 



353 . 



(3.58) 



where 61 > and 62 > are arbitrarily chosen under the constraint that ^3 > 0. 
Now we turn to the last three steps in the calculation : 



Step 1 : To deal with the first set of terms on the right hand side of (3.54) we use the expression 
for ||curl(^||2 in Lemma [s] with 5o = | and write 



Icurl C\^dV - (^\\v,,\\l + ^e^\\w,,\\i^ > |V2u;3p + \uyz\'^ + \vxz?]dV 



This turns (3.54) into 



< -2(53i?e-i / {|V2^^3|' + 
Jv 



(3.59) 



(3.60) 



'■yz\ 



+ \Vxz\'^]dV 



.-(1 - 2.3)i^e-3 /jV3e.Pdy + ( + ^) i.e3|H 



a 



1 + 



a 
64 



i?e-3||e||^ 



62 + 



1 

2&2 



\v\\i + ^m\l 



\l 



(3.61) 



Step 2 : Now this inequality is re-scaled back to dimensional variables defined in Table [2] and H^t) 
defined in (3.3). This involves multiplying both sides of (3.60) by Re^. 



,dH ^ 263L 
^ dt - 



+ 



^0 Jv 
9 243 



(I V2f23p + "al^il.xaxa P + aal^2,xiX3 H dV 



"hor II 6 



+ 



26f 12(5| 



(1 - 263)L^ 

— fi 2 

11^-11^ + -6^^'^ 64 



\^T,,^?dV 



1 



2^,-6 



0"a 



2b2) " 2,52 

Step 3 : Finally, we take the time integral over an interval [0, 

253L2 



(3.62) 




JV 



{ 



+ \V-^U{^))dVdt>\U{T^). (3.63) 
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Because of regularity [281 Ell 130', ^3J] ^.{t^,) is always under control from above and it also has a 
uniform lower bound 'H{t^,) > although zero may be a poor lower bound. 

Together with the use of Lemma [3] with = I. gives the result of Theorem [T] where Pu^or^ /^"s 
and Pi^T are defined in Table [5] and the forcing function Fh{Q) is defined in ( |3.2[ ). ■ 



4 Potential implications for simulations 

The main result of this paper is that solutions of HPE can potentially develop extremely small scales 
of motion, allowed by the estimates derived here. These size scales decrease as TZ~^ and 
which means they could easily become of the order of metres or less at the very large values of these 
parameters achieved in both atmospheric and oceanic flows. The hydrostatic estimate for the length 



scale defined \h in (3.7) is of the order of a metre or less. Of course, this very small estimate may 
not be the thickness of a front ; instead, it may refer to the smallest scale of features within a front. 
The importance of the tendency to produce vigorous intermittent small scales in NWP and ocean 
circulation simulations remains to be determined but it may effect parameterizations as numerical 
resolution improves. In particular, one may ask whether parameterizations developed at coarser scales 
will still be accurate at finer scales, if the finer scales undergo the extreme events whose potential 
appearance has been predicted in this paper. As for the perennial question of initial conditions, one 
must hope that flow activity initialized at coarse scales will be consistently followed to smaller scales 
without undue amplification of simulation errors. 

The consequences of Theorem 1 in ^is that space-time is potentially divided into two regions S+ 
and . The region could be a union of a large number of disjoint sets, and if it were non-empty 
the flows in would be dominated by strong concentrated structures. Very large lower bounds on 
double mixed derivatives of components of the velocity field {ui, U2, u-^) such as \d'^ui/dxsdx2\'^ 
or \d'^U2/dxidx3\^ may occur within thus representing intense accumulation in the (xi,X3)- 
and (x2, X3)-planes respectively. For a nonempty S~, one would see the spontaneously formation of 
front-like objects localised in space that would only exist for a finite time. TZuhor ^ /oca/ horizontal 
Reynolds number depending upon the local space-time values of u{xi,X2,X3,t) = y^nf -|- m| and 
Tla,T is a Rayleigh number dependent on the local temperature T(xi, X2, X3, r). The large lower 
bounds on double-derivatives of solutions within the regions can be converted into the large lower 
bounds on inverse length scales Xjj^. Thus to resolve a region such as this would requir^ 

Number of grid points > const {U^uhor + '^It + Re^'R-a,T) ■ (4.1) 

It is also worth remarking that the L^-norm arising in the proof of Theorem |3.ll leading to the sixth 
powers of the local Reynolds numbers Tlu^^r • precisely the norm that was proved by Cao 

and Titi [JSj [22] to be bounded for HPE. While TZu^or ^ function of space-time it is a bounded 
function, but how much Tluhor oscillates around its global space-time average Re is unknown : this 
could vary significantly in different parts of the flow. Thus, how = {ui,U2, u^} varies across a 
front is an important issue. The limitations of the result are that no further information is available 
from the analysis regarding the spatial or temporal statistics of the subsets of on which intense 
events would occur. 



72.U3 is expected to be negligible compared to TZu^^^ because U3 ~ 0{e) 
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If the regularity problem were to be settled in the NPE case, the results would likely be qualitatively 
the same but with a non-negligible TZu^ term whose contribution may be significant in regions 
of strong vertical convection. There would also have to be significant technical differences: the 
domain would need to made periodic in the velocity variables and their derivatives because of lack 
of specification of horizontal velocity derivatives. 

Future improvements in numerical capabilities for the prediction of weather, climate and ocean 
circulation may be expected to enhance spatial and temporal resolutions. In addition, they will raise 
the issue of the optimal allocation of numerical resources. For example, improving the computations 
for parameterizations of other currently unresolved physical processes (such as phase changes in 
cloud physics) may have effects that are at least as significant as computing non- hydrostatic effects 
at finer resolution. Improvements in resolution will also raise the issue of whether subgrid-scale 
parameterizations of these unresolved physical processes that have been developed for numerical 
prediction at coarser scales will transfer accurately to computations at finer scales, regardless of 
whether the hydrostatic approximation is retained. Thus, one may expect the HPE to remain central 
in the discussions about choices among the various potential numerical code implementations for 
weather, climate and ocean circulation predictions well into the foreseeable future. Even though they 
are mathematically well-posed, the HPE have been shown here to contain the potential for sudden, 
localized events to occur on extremely small scales in space and time. 
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